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Abstract 

In the 1920's, Madelung noticed that if the complex Schroedinger wavefunction is expressed in 
polar form, ip = i?e*^, then the quantities p = B? and u = V0 may be interpreted as the hydro- 
dynamic density and velocity, respectively, of a compressible fluid. In this paper, we generalize 
Madelung's transformation to the quaternionic Schroedinger equation. The non-abelian nature of 
the full SU (2) gauge group of this equation leads to a richer, more intricate set of fluid equations 
than those arising from complex quantum mechanics. We begin by describing the quaternionic 
version of Madelung's transformation, and identifying its "hydrodynamic" variables. In order to 
find Hamiltonian equations of motion for these, we first develop the canonical Poisson bracket and 
Hamiltonian for the quaternionic Schroedinger equation, and then apply Madelung's transforma- 
tion to derive non-canonical Poisson brackets yielding the desired equations of motion. These are 
a particularly natural set of equations for a non-abelian fluid, and differ from those obtained by 
Bistrovic et al. only by a global gauge transformation Because we have obtained these equa- 
tions by a transformation of the quaternionic Schroedinger equation, and because many techniques 
for simulating complex quantum mechanics generalize straightforwardly to the quaternionic case, 
our observation leads to simple algorithms for the computer simulation of non-abelian fluids. 
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I. INTRODUCTION 



Since the 1930's it has been known that it is possible to construct quantum mechanics 
over fields other than the complex numbers. In particular, a detailed and consistent presen- 
tation of quaternionic quantum mechanics has been developed and found to have a number 
of fascinating theoretical properties 0. In the present work, we develop a relationship be- 
tween quaternionic quantum theory and the simulation of non-relativistic, non-abelian fiuid 
dynamics. Many computer simulation techniques developed for complex quantum mechanics 
generalize straightforwardly to the quaternionic case, so the observation made in this paper 
makes possible new simulation strategies for non-abelian fiuid dynamics. 

The Madelung transformation transforms the complex Schroedinger equation into the 
continuity and Euler equations describing an inviscid fiuid. The magnitude squared of the 
wavefunction is the mass density, and the gradient of the phase of the wavefunction is the 
hydrodynamic velocity. The mass continuity equation follows directly from the conservation 
of probability density, and the momentum equation of motion follows with the identification 
of a density dependent pressure. 

Recently, there has been interest in developing non-abelian fiuid dynamical models to 
describe systems such as quark-gluon plasmas, where one may wish the color degrees of 
freedom to be retained in the fiuid description [JI], ^ |^. The color degrees of freedom are 
associated with the gauge symmetry group of the field; for the example of a quark-gluon 
plasma this would be the group SU{3) of quantum chromodynamics (QCD). 

A natural question to ask in the context of quaternionic quantum mechanics is: What is 
the generalization of the Madelung transformation for the quaternionic Schroedinger equa- 
tion? As the quaternionic elements are anticommuting, and in particular have the Lie algebra 
of SU(2), we might plausibly expect the fiuid equations arising from such a transformation 
to be those of a non-abelian fiuid. 

The remainder of this paper is organized as follows. We first identify the quaternionic 
generalization of the Madelung transformation, and identify the hydrodynamic variables. In 
order to find evolution equations for these quantities we cast the quaternionic Schroedinger 
equation in Hamiltonian form, defining a canonical Poisson bracket and Hamiltonian for 
the quaternionic wavefunction. We then use our identification of the quaternionic Madelung 
variables to define a non-canonical transformation of our Hamiltonian field theory, obtaining 



2 



the bracket structure and evolution equations. We compare our approach with that of Jackiw 
et al. 0, ^ , and briefly present several numerical algorithms by which our equations may 
be straightforwardly simulated. We close the paper with some conclusions and directions 
for future work. 

II. MADELUNG VARIABLES 

The identification of the gradient of the phase as the relevant hydrodynamic variable in 
the complex case is a strong indication that local gauge symmetry plays an important role in 
the transformation, and will therefore be useful as a guide in the more difficult quaternionic 
case. We first identify the gradient of the phase of the wavefunction, which is identified 
with the fiuid velocity in the complex Madelung transformation, as a pure U{1) gauge. The 
condition that the fiuid velocity be a pure gauge translates to the constraint of vanishing 
vorticity for the Madelung fiuid. We then use this observation to obtain the form of a pure 
gauge for a quaternionic SU (2) Yang- Mills gauge field. We find three vector fields and 
show that they obey constraints which are the generalization of vanishing vorticity to the 
non-abelian case. 

To perform the program outlined above we first construct the gauge transformation of a 
gauge field Ai with arbitrary gauge symmetry group Q. Our gauge field is introduced via 
the gauge covariant derivative Da = da — Aa- All physical quantities are invariant under 
the transformations 

^ ^ ^' = # (1) 

and 

A'^, (2) 

where g is an element of the gauge symmetry group and the transformation properties of Aj 
are obtained by requiring 

{da-A'^)gij = g{da~Aa)iJ, (3) 

which gives in turn 

A'^={dag)g-' + gAag-\ (4) 

For complex quantum mechanics coupled to the electromagnetic field the gauge group 
Q is U{1), the element g is any unimodular complex number, and the quantity {d^g^g~^ is 
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simply the gradient of the change in phase of the wavefunction. The above transformation 
must produce no change in electromagnetic field strength Fap = daAp — dpAa- It is clear 
that they do not, as the change in takes the form of a similarity transformation and the 
addition of the pure gauge {dag)g^^- The condition that the gauge fields remain unchanged 
is the elementary condition V x V^^ = 0, i.e., the constraint that the Madelung fluid has 
vanishing vorticity. 

The above discussion motivates a straightforward generalization of the Madelung transfor- 
mation. The Madelung velocity fields have the form of pure f/(l) gauges, and the condition 
that pure gauges do not contribute to the gauge field strengths translates into the condition 
of vanishing vorticity for the Madelung velocity field. Turning to the quaternionic case, 
the gauge group Q is SU{2), the element g is any unimodular quaternion, and we identify 
the pure gauges {dag)g~^ with the relevant hydrodynamic variables for the quaternionic 
Madelung transformation. 

We write the group element g as 

g = e'^'e^^^" (5) 

Where z, j and k are the quaternionic elements satisfying i"^ = = k"^ = ijk = —1, and /i, 
u and 6 are real angles. Substitution of this form for g into Eq. gives 

( V^) g-^ = iVfi + e^f'je-'^'Ve + e^f^e^^ke'^^e-'^'Vu, (6) 

We write 

(V^)^"' = ^u + jv + A;w, (7) 

where 

u = Vyu + COS 2eViy 

V = cos2/iV6' + sin26'sin2/iVi/ (8) 
w = sin — sin 26 cos 2/i Vi^. 

The definition of the density is p = ipip, where ip is the quaternionic wavefunction ip = 
y/pg and its quaternionic conjugate ip is ip = \fpg~^ = y/pe'^^e'^^e'^^. The quaternionic 
generalizations of the Madelung variables are then p, u, v and w. 

These quaternionic Madelung variables are the pure gauges of an SU{2) Yang-Mills the- 
ory. We can see this by showing that the Yang-Mills field strengths are zero. The non-abelian 
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generalization of the Maxwell field strength is 

F^p = d^Ap - dpA^ - Ap] (9) 

where A^ are quaternionic gauge fields. Sufficient and necessary conditions for these field 
strengths to be zero are 

V X u = 2v X w 

V X V = 2w X u (10) 

V X w = 2u X V 

It may be readily verified that these equations are satisfied by u, v and w, and hence are the 
quaternionic generalization of vanishing vorticity in the complex Madelung transformation. 
We now obtain the equations of motion for u, v and w. 



III. EQUATIONS OF MOTION 

In order to obtain the equations of motion for u, v and w we first obtain a canonical 
Hamiltonian form for the quaternionic Schroedinger equation. To do this we need to define 
a Hamiltonian and a Poisson bracket. Defining the Hamiltonian is straightforward, as it 
follows directly from the definition for the complex case, 

H = - I i)Hi)d^z = - I ( --V^VV + Vijij)(fz 

^ 2 ^ (11) 

1 /"/^IVPP , P^^2^^2^^2^^y\ ^3^^ 



2 J \ 8p 2 

where we use units such that h = m = 1 throughout. By the argument presented in 
Appendix ^ we define the bracket 



5Q(x)5P(y) 5g(x)5P(y)\ 3 



^ ' ' 6d{z) 6c{z) 6c{z) 6d{z) '"^ ^' 
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where a, b , c and d are the quaternionic components of the wavefunction ip = a + ib+jc+kd. 
The above Poisson bracket and Hamiltonian reproduce the quaternionic Schroedinger equa- 
tion. The definitions Eq. (P) are then a non-canonical transformation of this Hamilto- 
nian field theory. The bracket structure of the quaternionic Madelung equations follows by 



lengthy but straightforward algebra. We spare the reader the details and present the non- 
canonical brackets amongst p, u, v and w in Appendix ^. Using the above Hamiltonian, 
the equation of motion for p is then 

dp 



dt 



+ V ■ (pu) = 0. 



(13) 



This is the continuity equation, and if p is identified as the singlet density, this equation 

identifies u as the singlet velocity. The equation of motion for u is 

du , 1 

— = -W - -' 

dt 2 

where V' is the density-dependent potential (also called the Bohmian potential) which also 
appears in the complex Madelung transformation. 



w 



V w 

-V-(pv)--V-(pw), 



(14) 



v = v - 



:V 



Vp 



(15) 



Note that the curl of Eq. (0) is generally nonvanishing. Since u is identified as the singlet 
density, this means that the quaternionic Madelung fluid has non-zero vorticity. This can be 
clearly seen by inspecting the {u(x), u(y)} bracket, which has the correct form for rotational 
flow. The equations for v and w are 



and 



dt 2 



9w 
~dt 



'V ■ (pw)" 


- 2w 


P 





rv-(pv)i 


+ 2v 


. p . 





y'(x) + ^(u2 + v2 + w2) 



V'(x) + ^(u2 + v2 + w2) 



u. 



+ -V ■ (pw) 



u 



p 



V-(pv). 



(16) 
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One immediately notices that the equations for v and w are invariant under the duality 
transformation 

v' \ / cos Q — sin ^ \ / v 
w' / \ sin Q cos 9 I \w 
This duality symmetry is a consequence of the global gauge invariance of the quaternionic 
Schroedinger equation. This global invariance is discussed in detail in Appendix y. 
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IV. NON-ABELIAN FLUID DYNAMICS 



Recently there has been interest in constructing non-abelian fluid dynamics and fluid mag- 
netohydrodynamics in order to simulate the long-wavelength degrees of freedom in a quark- 
gluon plasma. In particular, Jackiw et al. have described several models for the relativistic 



dynamics of such a system p, Q . Those authors have also presented a non-relativistic multi- 
component wavefunction very similar to the quaternionic wavefunction considered here [|1|. 
In fact, as we shall show, their model differs from the quaternionic Madelung transformation 
we have investigated only by a global gauge of the type described in detail in Appendix y. 
In this section we compare their multicomponent wave equation with our approach and give 
the global gauge transformation relating their work to our own. 
Bistrovic et al. consider a two component complex wave equation, 

.^^4V=.. ,19, 

Those authors consider the group SU{2) and the fundamental representation: = a"'/ (2i), 
{T", T^} = — (5"^/2, where cr", a = 1,2,3 are the Pauli matrices. The T"'s here are simply a 
two dimensional complex representation of the quaternionic elements 

Ti = |,r2 = |,r3 = |. (20) 

Those authors define the conserved singlet and color current; 

f = {cp,i), P = V'V, j = ImV'^V^. (21) 

= {Cpa, J a) , Pa = #^T> , Ja = Re ^^T'^V^J . (22) 

and the wavefunction, 

^ = VpgA (23) 

where p is the scalar tp'^ip, where f denotes hermitian conjugation, g is a. group element, and 
A is an arbitrary constant vector taken to be 

(24) 

in which case iA'''T"A = 5°'^/2. The singlet density is p, while the singlet current j is 

j=pp, p = -tA^ g-'Vg A. (25) 

The key relation connecting our approach to Jackiw's is: 

g-'Vg = -2p^T'' (26) 
p = p"t", (27) 




where = iu^T°-u = (5"^/2, and hence 



p = p'^t" = p'^. (28) 

It follows that the singlet velocity in is equated to the negative of third quaternionic com- 
ponent of the wavefunction in our approach. His multicomponent wavefunction is therefore 
equivalent to a quaternionic Schroedinger equation 

This quaternionic Schroedinger equation is related to the one considered in this paper by a 
one-parameter family of global gauge transformations, described in detail in Appendix p. 

The color and singlet currents identified in [|l| may be obtained from the fields u, v, and w 
defined for our quaternionic Madelung equation. The quaternionic Madelung transformation 
we have described is a particularly simple implementation of this type of non-abelian fluid 
model, but its usefulness is limited as the approach cannot be extended to groups other 
than SU{2). The quaternionic approach has some advantages in terms of ease of numerical 
simulation, which we describe below. 

V. SIMULATION 

The Schroedinger equation may be simulated by a variety of methods. We briefly describe 
three methods, all of which may be trivially generalized to the simulation of quaternionic 
quantum mechanics. The first method is operator splitting, or its higher order variant, 
the symplectic integrator method. This method is perfectly unitary, can be generalized to 
three dimensions, and can be generalized to high orders of accuracy |^, 0, H, ^, |T0|. The 



second-order spatial derivatives required in the operator splitting method may be avoided 
by using a quantum lattice-Boltzmann method (QLBE) |Tl|, |I2|, [1^. The QLBE is uncondi- 
tionally stable, perfectly unitary, and offers considerable advantages for parallel computing 
applications. 

Closely related to the quantum lattice Boltzmann approach is the quantum lattice-gas 
model. These algorithms were originally developed for simulation of the Dirac equation 



in one spatial dimension |15|, |16|, [T^- The non-relativistic quantum lattice gas was 



subsequently developed for the simulation of many-body quantum mechanics in arbitrary 
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numbers of dimensions [0. The quantum lattice-gas algorithm for a single nonrelativistic 
particle provides a tractable technique for the simulation of the quaternionic Schroedinger 
equation. 

It should be noted that the quantum lattice gas algorithm (QLGA) of was developed 
in order to find a quantum algorithm for Schroedinger's equation. The QLGA provides 
a method of simulating quantum mechanics on a quantum computer with an exponential 
enhancement in the performance. Minor modifications to this QLGA should also provide a 
method of simulating the quaternionic Madelung equations derived above on a (complex) 
quantum computer with a similar exponential speedup. 

VI. CONCLUSIONS 

We have obtained the quaternionic generalization of the Madelung transformation. We 
have identified the generalizations of the hydrodynamic fields and obtained their evolution 
equations by a non-canonical transformation of the Hamiltonian field theory for the quater- 
nionic Schroedinger equation. 

The equations we have obtained also possess a limit which may be of interest. The 
choice of arbitrarily large potential may always make the second terms in Eqs. (|l^) and (|l^ 
dominant. In this case a solution which is constant in space and oscillatory in time will 
be a good first-order approximation to the behavior of v and w. This suggests that there 
may exist an approximate treatment of these equations in which v and w are treated as 
"fast" variables which may be averaged away. This would be analogous to the gyro-averaged 
description of MHD, in which the fast gyrofrequency motions are averaged away and only 
the slow motion of the gyro centers is retained in the reduced dynamics |I9| , |20| , pT| . The 



Hamiltonian approach taken in this paper provides a good starting point for the derivation 
of such an averaged treatment. 

The Madelung transformation of complex quantum mechanics seems intially computa- 
tionally promising, as it shows us how to simulate the nonlinear Euler equations by sim- 
ulating the linear Schroedinger equation. However, irrotational inviscid fluid flow is not a 
problem of sufficient interest to merit serious computational investigation. The quaternionic 
Madelung transformation gives rise to equations of considerably more interest, as they de- 
scribe the (rotational) flow of a non-abelian fluid. However, the constraints in Eq. ([T0|) are 
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the restriction on the non-abehan fluid, analogous to vanishing vorticity in the abehan case. 
This suggests that there exists a generahzation of the equations obtained here to non-abehan 
flow flelds which do not satisfy the constraints of Eq. (p^Ol). It is possible that these general 
non-abelian flow equations will contain as much additional physics as does rotational Euler 
flow compared with irrotational Euler flow. 

A natural question arises as to whether one may generalize this approach to groups other 
than SU (2) and flelds other than the quaternions. We may answer this immediately in the 
negative: There are no division algebras other than the reals, the complex numbers, the 
quaternions and the the octonions [^. The octonions form a division algebra, but are not 



associative, and this precludes the formulation of an octonionic quantum mechanics 0. 

Finally, we have noted that the computational tractability of the complex Schroedinger 
equation is shared by its quaternionic generalization, and so we anticipate that extensive 
simulation of the quaternionic Schroedinger equation by any one of the methods described 
above will yield results of interest to the particle physics community. 
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APPENDIX A: QUATERNIONIC POISSON BRACKET 

If we write the Schroedinger equation as a pair of real equations for the real and imaginary 
parts of the wavefunction ip = a + ib v/e obtain: 

(Al) 

where H = — + V. Note that the above construction does not enable us to construct 
quantum mechanics over the reals. The symplectic matrix in the above equation does not 
have real eigenvalues and so there are no energy eigenstates in real quantum mechanics pf. 
However, the form of this equation makes it immediately obvious that the Schroedinger 
equation may be described by a canonical Hamiltonian field theory in which the real and 
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imaginary parts of the equation are canonically conjugate. The definition of the bracket 
follows automatically, 

'(5g(x)5P(x) 5g(x)5P(x) 



{Q(x),P(y)} 



(A2) 



6a{z) 6b{z) 6b{z) 6a{z) 

The quaternionics have a four dimensional real representation, which enables us to write the 
quaternionic Schroedinger equation as four real equations for the four components of the 
wavefunction ip = a + ib + jc + kd, 

/ ^0 \ 

-H 
-H 
\ H J 

Again, we see that these equations have a symplectic structure where a and b are canoni- 
cally conjugate, and so are d and c. The definition of the bracket, Eq. (|T2]), again follows 
automatically. 









d 


b 




~dt 


c 






\d) 





b 
c 
\d/ 



(A3) 



APPENDIX B: POISSON BRACKET STRUCTURE 

Our quaternionic Madelung variables have the following bracket structure: 

{p(x), Wi{y)} = -4(5(x - y)vi{yi) 

{p(x),t;,(y)}=45(x-yH(x) (Bl) 
{p(x),t.,(y)} = -2(5:(x-y) 
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{«i(x),Mj(y)} 

{ui{x),Vj{y)} 

{vi{x),Uj{y)} 
{vi{x),Vj{y)} 



{vi{x),Wj{y)} 
{wi{x),uj{y)} 
{wi{x),Vj{y)} 
{wi{x),Wj{y)} 



2(^(x-y) 



25^(x-y)i;,(x)^45(x-y) 



P(x) 



P(x) 



-(aiVj(x) - ajVi(x)) - Mi(x)Wj(x) 



25'(x-y)w,(x) 4<5(x-y) 



P(x) 



P(x) 



1 



(5jWi(x) - i9iWj(x)) - i;j(x)Mj(x) 



25Ky-xK-(y) 45(y-x) 



p(y) 



p(y) 



25,'(y-x)M,(y) , 25j(x-y)M,(x) 



p(y) 



+ 



P(x) 



45(y - x) 



P(x) 

5^;.(x-y) 5^(x-y)9,p(x) 



- ( a.K,(x) - a,u.(x) ) + ..(x)^,(x) - .,(x)^.(x) 



P(x) 



P(x)' 



45(y - X 



Mj(x)^j(x) 
P(x) 



2<5Ky-xK(y) ^ 4c5(x-y) 

p(y) p(x) 
^.'^•(y-^) , <^^(y-x)9,p(y) 



-{diWj{x) - djWi{x.)) - Vi{yi)Uj{x) 



p(y) 

5.'(y-x) 



+ 



= -2 



p(y) 



+ 45(x-y) 
Wi(x) 



^^(x)^j(x) 
P(x) 



+ 45(y-x: 



p(y) 



;(5i%(y) - 9j^ii(y)) + Vi{y)wj{y) - Vj{y)wi{y) 



(B2) 



APPENDIX C: GLOBAL GAUGE SYMMETRY 



We obtained our quaternionic Madelung variables by considering the local gauge invari- 
ance properties of our Schroedinger equation coupled to a quaternionic gauge field. The 
observables of quaternionic quantum mechanics posess global SU (2) gauge invariance, and 
this invariance will lead to a set of symmetry properties for the quaternionic Madelung 
equations. The Schroedinger equation is 



Under the global gauge transformation ip ^ ip' — (pip this becomes: 



(CI) 



(C2) 
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If is a unimodular quaternion then the transformation (f)~^Q(j), where Q is any quater- 
nion, is an automorphism of the quaternions, and (j)~^Q(f) is also a unimodular quaternion. 
The equivalent global gauge transformation applied to the complex Schroedinger equation 
also induces an automorphism of the complex numbers, but in that case it is the trivial 
automorphism V^^ : e^*^Ze*^ = Z. The choice of i in Eq. (PT|), which we inherited from 
the complex Schroedinger equation, constitutes a global gauge choice. We could equally 
well write the quaternionic Schroedinger equation replacing i by j, or fc, or indeed any pure 
imaginary unimodular quaternion. We now consider how such global gauge choices affect 
our quaternionic Madelung equations. 
Our Madelung variables are defined by: 

{Vg)g-^ = iVL + jw + k^, (C3) 

applying a global gauge transformation to this quantity gives: 

0"^(V^) g-^(t) = (j)-H(f)u + 0"^j0v + (j)-^k(j)w, (C4) 

Taking the quaternion automorphism defined by: 

fp = (pCpcf). (C5) 

Where Cp = i,j,k,p = 1,2,3 are our original quaternionic elements and fq are our new 
quaternionic elements. We express the new quaternionic elements in terms of the old: 

We can determine the properties of the coefficients a^p from the algebra obeyed by both sets 
of quaternionic elements: 

CqCrf Sq^ ~1~ ^qrs^s (^'^) 

Substituting: 

Oimp^nq I ^pq ~\~ ^pqs^^s 1 ^mn ~\~ ^mno^or^^r (^^) 



Where e^no is the completely antisjnumetric Levi-Civita tensor. Separating the coefficients 
of the quaternionic elements gives: 



Cirp^sq^pq 



_ (C9) 
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These equations show that the matrix apq is orthogonal with unit determinant. The global 
gauge transformation induces a linear transformation (actually a proper rotation) of our 
Madelung variables, 



u = Oiiu + a2iv + flsiw 
v' = ai2U + a22V + a32W 
w' = ai3U + a23V + asaw. 



(CIO) 



The Madelung transformation therefore applies to a three-parameter family of quaternionic 
Schroedinger equations, and the results of simulating any member of this family may be 
related to the results of simulating any other member by a global gauge transformation. 
As a specific example of this we construct the transformation relating the quaternionic 
Schroedinger equations: 



and 



.(hp_ 

''Itt 



dt 



(Cll) 



(C12) 



The two equations are related by the global gauge transformation ip' 



where 



Ad 



2, and the matrix apq is 



"pq 



(C13) 



-A 

0-10 
-10 

This is the global gauge transformation relating the quaternionic Schroedinger equation 
considered here to the multicomponent non-abelian Schroedinger equation considered by 
Bistrovic et al. 

The gauge transformations relating the two quaternionic Schroedinger equations above 
form a one parameter family. There is always a one-parameter subgroup of the full gauge 
group which leaves any quaternionic Schroediger equation identically invariant. For the 
Schroedinger equation, Eq. ( |(J11| ), the elements of this group may be written g = e*^/^. This 
subgroup induces the transformations 

1 ^ 

cos 9 — sin e (C14) 
sin 6* cos 6* 



"PQ 



I 
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This is the set of global gauge transformations which give rise to the duality symmetry of v 
and w noted in Eq. (p^). 
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